Spinons and helimagnons in the frustrated Heisenberg chain 
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We investigate the dynamical spin structure factor S{q,(jo) for the Heisenberg chain with ferro- 
magnetic nearest (Ji < 0) and antiferromagnetic next-nearest (J2 > 0) neighbor exchange using 
bosonization and a time-dependent density-matrix renormahzation group algorithm. For |Ji| <C J2 
and low energies we analytically find and numerically confirm two spinon branches with different 
velocities and different spectral weights. Following the evolution of S{q,oj) with decreasing J1/J2 
we find that helimagnons develop at high energies just before entering the ferromagnetic phase. 
Furthermore, we show that a recent interpretation of neutron scattering data for LiCuV04 in terms 
of two weakly coupled antiferromagnetic chains (| Ji| ^ J2) is not viable. We demonstrate that the 
data are instead fully consistent with a dominant ferromagnetic coupling, J1/J2 ~ —2. 
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Frustrated spin models often show complicated phase 
diagrams. Apart from phases with conventional (quasi) 
long-range magnetic order, also valence bond solid, chi- 
rally ordered, multipolar as well as spin liquid phases 
are possible. In addition to identifying and studying the 
different phases, also the quantum phase transitions be- 
tween them have attracted considerable interest [1, 2]. 
Numerical studies are hindered by the so-called sign 
problem making quantum Monte- Carlo simulations im- 
practicable. Our knowledge about frustrated systems in 
dimensions d > 2is therefore still limited [3-5] . In one di- 
mension, on the other hand, bosonization techniques can 
be used to analytically study the weakly frustrated case 
[6, 7]. Furthermore, with the density- matrix renormal- 
ization group (DMRG) [8] a powerful numerical method 
is available to calculate static properties at zero [6] and 
finite temperatures [9] . Lately, it has been shown that it 
is a useful strategy to study frustration in two or three 
dimensions by a successive coupling of one- dimensional 
chains [3, 10] thus further motivating investigations of the 
static and dynamic properties of frustrated spin chains. 

The standard spin model to study frustration in one 
dimension is the Ji — J2 Heisenberg model 

H = ^ [JiSjSj^i + J2SjSj^2] (1) 
j 

with antiferromagnetic coupling J2 > 0. Here Sj is a 
spin- 1/2 operator acting at site j. This model has re- 
cently been studied intensely for ferromagnetic coupling 
Ji < [11-15], driven by the discovery of edge-sharing 
cuprate chain compounds for which Eq. (1) seems to be 
the minimal model. These cuprates show fascinating 
properties including multiferroicity [16-20], i.e., an in- 
tricate interplay between incommensurate magnetic and 
ferroelectric order [21]. It is, however, so far unclear if all 
of their magnetic properties can be reasonably well de- 
scribed within the simple minimal model (1). In one of 
the best studied edge-sharing chain cuprates LiCuV04, 



for example, susceptibility data seem to point to a dom- 
inant ferromagnetic coupling a = J1/J2 ~ —2 [9, 22] 
while neutron scattering data have been interpreted in 
terms of two weakly coupled antiferromagnetic chains 
a ^ —0.7 [23], a conclusion which has later been chal- 
lenged [24-26]. 

Apart from being relevant for the edge-sharing cuprate 
chains, the dynamical properties of the Ji — J2 chain are 
also of fundamental interest. For a = the model con- 
sists of two decoupled antiferromagnetic chains whose el- 
ementary gapless excitations are spinons. Introducing a 
small coupling Ji < between the chains leads classically 
to a spiral magnetic order while bosonization predicts 
incommensurate spin correlations and an exponentially 
small gap in the quantum S = 1/2 case [7]. At ac = —4 
there is an unusual quantum critical point [2] separating 
the incommensurate from a ferromagnetic phase. Chang- 
ing the frustration ratio a < thus turns antiferromag- 
netic spinons through an incommensurate phase into fer- 
romagnetic magnons. 

In this letter we present a systematic study of the dy- 
namical spin structure factor 

^(^'^) = ^ E^"''^'"'"^ / dte'-\S]{0)S].{t)) (2) 

of the J1-J2 Heisenberg model (1) from the limit = 
of decoupled antiferromagnetic Heisenberg chains, across 
the quantum critical point o^c, into the ferromagnetic 
phase, a < —4. Note that due to SU{2) symmetry it 
is sufficient to consider the longitudinal correlation func- 
tion in (2). For \a\ <C 1 we compare our data with results 
obtained by bosonization while for ^ —4 we compare 
with spin wave theory. Finally, we present a comparison 
of our results with recent neutron scattering data for the 
multiferroic cuprate LiCuV04 [23]. 

We start by considering the weak coupling limit |a| <C 
1 by bosonization. On each of the two antiferromagneti- 
cally coupled sublattices we write 
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Here ^1(2) are bosonic fields obeying ttie standard com- 
mutation rules [(j)a{x)^dxOa'{x')] = i5a,a'^{x — x') and 
K12 are the Luttinger parameters. Let us first consider 
the free fermion case with J2SiSj ^{Sf Sj + h.c). 
Ignoring irrelevant terms, bosonization leads to 



(4) 



for each of the chains a = 1, 2. In this case Ki = K2 = 1 
and vi^2 = with vp = 2J2 in units of the lattice 
constant due to a doubling of the unit cell. Apart from 
irrelevant terms, the interchain coupling Ji introduces a 
density- density type interaction 
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We can absorb this term into the Gaussian part (4) by 
defining the new fields = (01 ±02)/ V2 and 0± = {Oi± 
^2) /a/2- The Hamiltonian (4) stays invariant under this 
transformation with a = 1, 2 being replaced by = +, — 
but the velocities are renormalized with vi^2 = 
vf ^ Ji/tt. 

In the isotropic Heisenberg case for a = the low- 
energy properties are still described by (4) but with 
vi^2 = ^27r and K12 = 1/2 as known from Bethe ansatz. 
To keep the SU{2) symmetry intact, we now use non- 
Abelian instead of the Abelian bosonization (3) for the 
interchain interaction Ji. This leads to various marginal 
terms [7, 27] and, in particular, a term which can again 
be expressed as in (5) but with a different prefactor 
2Ji/7r^3Ji/7r leading to a velocity renormalization 



vi^2 = 7TJ2 v± = 7rJ2(l lb 3ce/27r^) . 



(6) 



The other marginal terms induced by the coupling Ji 
have been studied by renormalization group methods 
in [7] and seem to lead to an exponentially small gap. 
Numerically, however, no gap in the parameter regime 
— 4 < a < has been confirmed yet and we will therefore 
neglect these terms for the moment. 

At weak coupling Ji, the structure factor S{q,iu) will 
have low-energy contributions at q ~ 0,7r/2,7r. We now 
calculate these contributions in the isotropic case using 
the free boson Hamiltonian (4) with = +, — and ve- 
locities v± given by Eq. (6). Using the identity (3) — 
rewritten in terms of (t)± = (01 zb 02)/\/2 — and the stan- 
dard result for the boson propagator we obtain 

S{q>0[7^-q<7r],Lo) (7) 
^ Kq{l zb cos q)S{(jU — Vj^q) + Kq{l =F cos q)6{(jj — v^q) 

q^ 

^ 2KqS{uj - v^l_]q) ^ K—S{uj - v_i^]q) . 



The singularity ai u = V-^q dominates S{q ~ O^u) while 
the singularity at uj = v^q is dominant for g ^ tt. The 
other contribution is in both cases surpressed by a factor 
g^/4. Note that V-^ < V- because Ji < 0. Irrelevant 
band curvature and interaction terms will lead to a finite 
linewidth [28, 29], an aspect which is beyond the scope 
of the present study. For \k\ = |7r/2 — g| <C 1 and 
sufficiently different we find 
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Here K± ^ 1/2 are the Luttinger parameters for the 
two modes at weak coupling Ji . Note that the divergen- 
cies in (8) are weaker than the square-root divergence 
for the pure isotropic Heisenberg chain [30]. Taking the 
additional marginal current- current interactions into ac- 
count perturbatively, we find furthermore an asymmetry 
S{7r/2 - q,uj) > S{ti/2 + q, uj) with g > smah. 

To test the analytic predictions in the weak coupling 
limit and extend our study to the experimentally relevant 
regime of strong interchain coupling, we now turn to a 
numerical calculation of the dynamical spin structure fac- 
tor. We use an adaptive time-dependent density-matrix 
renormalization group (DMRG) algorithm with a second 
order Trotter- Suzuki decomposition of the time evolution 
operator [31]. As time step we choose J2(5t = 0.1. We 
present results for an open chain with N = 400 sites 
with 400 states kept in the adaptive Hilbert space. In 
order to perform the Fourier transform in time in Eq. (2) 
one has to deal with the problem that numerical data 
are only available for a finite time interval t G [0,tmax]- 
The maximal simulation time, tmax ^ 40 J2, up to which 
our numerical results are reliable has been estimated by 
keeping track of the discarded weight and by comparing 
with exact results for the XX model and Bethe ansatz 
results for the isotropic Heisenberg chain [32]. We calcu- 
late the spin correlations {S^^2i^)'^N/2±j(^))- Since two- 
point spin correlations are negligible for distances much 
larger than v±t^ we obtain results almost unaffected by 
the boundaries in the accessible time interval. These data 
are then extended in time using linear prediction [33, 34] 
leading to a smooth exponentially decaying extrapolation 
of the data. We want to stress that linear prediction does 
not allow to obtain reliable results for t > tmax but rather 
represents a smooth cutoff which does not affect the data 

for t < tmax- 

In Fig. 1 we show our results for S{q^uj) for various 
frustrations. A further check of the quality of the numer- 
ical data is obtained by considering the sum rules Ii = 
= (5,^5^,) where 5| = Ar-V2^.e-«5| 
and I2 = ^qli = 1/4. For all frustrations shown we 
find that the sum rules are fulfilled with an absolute er- 
ror of at most 2%. As an example, we show results for 
/i at a = -2 in Fig. 2(a). 

For weak frustrations a = —1/2,-1 (Fig. l(a,b)) two 
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FIG. 1: (Color online) S{q,(jj) for frustrations (a) a = —1/2, (b) a = —1, (c) a = —2, (d) a = —3.5, (e) at the quantum critical 
point a = —4, and (f) in the ferromagnetic regime a = — 5. The insets to (b) show a constant q and constant cj scan with 
arrows denoting spinon modes discussed in the text. The circles/dashed lines in (d-f) denote the magnon dispersions, Eq. (9). 



excitations with different velocities v± are clearly visi- 
ble near q it/2. As expected from bosonization, the 
smaller velocity v+ agrees with that of the dominant ex- 
citation at ~ while the one at g ~ tt has velocity 
V-. The velocities extracted from the numerical data 
also agree fairly well with the prediction from bosoniza- 
tion, see Fig. 2(b). A closer inspection of the data for 
a = —1 also reveals the second mode with velocity V- at 
q ^ {) with a weight suppressed by approximately q^ /A 
relative to the dominant mode, see inset of Fig. 1(b). 

At larger frustrations the weak coupling picture from 
bosonization clearly breaks down. For a = —2, —3.5 
shown in Fig. 1(c) and (d) the low energy spectral weight 
is concentrated at an incommensurate wave vector qi. 
Classically, frustration leads to the formation of a he- 
lical state with a pitch vector Q = arccos(|a|/4). For 
the quantum model, incommensurate spin correlations 
have been shown to occur with wave vectors qi which 
approach 7r/2 with increasing a much faster than in the 
classical case [9, 11] in full agreement with our dynami- 
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FIG. 2: (Color online) (a) Sum rule Ii for a = — 2 from 
frequency-integrating the numerical data (symbols) compared 
to a static DMRG calculation (line). The inset shows the 
absolute error of the frequency- integrated data, (b) Velocities 
of the elementary excitations v±. The symbols denote the 
values extracted from the numerical data for S{q,u) at q ^ 
0,7r the lines the prediction from bosonization, Eq. (6). 



cal data. For a = —3.5 we observe — in addition to the 
low-energy spectral weight near qi ^ 0.267r — the develop- 
ment of three magnon-like dispersions at higher energies. 
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For the classically expected state with long-range spiral 
order at wave vector Q, spin- waves have the dispersion 

6, = ^Al-Bl (9) 

with Aq = -Jq + JJ2 + (Jg+Q + Jq-q)/^. Bq = 

-Jg/2 + (Jg+Q + Jg_Q)/4 where Jq = Ji cos (7 + J2 cos 2(7 
[35, 36]. In the quantum model, long-range order 
is destroyed. However, close to the quantum critical 
point ac = —4 where the model starts to develop a 
long-ranged ordered ferromagnetic state, the correlation 
length will become large so that the excitations at high 
energy remain helimagnon-like. As shown in Fig. 1(d) 
the helimagnon dispersion (9) does indeed describe the 
high-energy modes very well with a renormalized effec- 
tive a = —3.85 corresponding to a pitch vector Q = 
arccos(|a|/4) ^ 0.27. The three modes [35] are then 
given by Cq and eq±q. where Qi is the incommensurate 
wave vector of the quantum model. 

For a < — 4 the ground state is a simple ferromag- 
net, (5 = 0, and (9) reduces to the magnon disper- 
sion Cq = —J2 [(1 — cos 2q) + a{l — cos q)] which is in ex- 
cellent agreement with the numerical data as shown in 
Fig. 1(f). At the quantum critical point, ac = —4, the 
magnon dispersion becomes quartic, Cq = J2^^/2 at small 
q. As already noticed in [2] spin- wave theory does not de- 
scribe the quantum critical point correctly at low energies 
due to the degeneracy of the ferromagnetic with valence 
bond solid states [37] . This is confirmed by our data (see 
inset of Fig. 1(e)) showing that the low-energy spectral 
weight does not follow the magnon dispersion while the 
agreement is good at higher energies. 

Let us finally discuss our results in the light of re- 
cent neutron scattering experiments on LiCuV04 [23] 
and LiCuSb04 [38]. The experimental data for LiCuV04 
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FIG. 3: (Color online) Neutron scattering data taken from 
Fig. 3 of Ref. [23] (circles) compared to the DMRG data with 
J2 = 6 meV and a = — 2 (solid lines). The dashed lines 
are the same DMRG spectra broadened by convoluting with 
a Gaussian with a FWHM of 1.1 meV. The height of the 
theoretical spectra has been scaled to fit the data best. 

have been interpreted in terms of the J1-J2 Heisenberg 
model with |Ji| < J2. However, none of the features 
predicted by bosonization in this limit and visible in 
Fig. l(a,b) have been observed. The neutron data in 



Fig. 2 of [23] look instead remarkably similar to our nu- 
merical results for a = —2. Comparing to the numeri- 
cally calculated structure factor one has to keep in mind 
that the neutron scattering data are slightly supressed at 
higher momenta due to the g'-dependence of the atomic 
form factor for Cu^+ ions and the whole spectrum is 
broadened due to a finite q, uo resolution. However, this 
does not affect the qualitative features of the spectrum. 
In fact, we can even obtain a quantitatively satisfying 
description of the data, see Fig. 3. For small momenta 
the agreement becomes excellent when convoluting the 
numerical data with a Gaussian to take the finite instru- 
mental resolution into account. A frustration a ~ — 2 
implies an incommensurate wave vector qi ^ 0.457r, see 
Fig. 2(a), consistent with the neutron data and suscepti- 
bility measurements [9, 22]. For such a strong frustration 
the explanation of the spectral weight at high energies 
in terms of multispinon excitations of the antiferromag- 
netic chains offered in Ref. [23] is not viable. Instead, this 
weight is related to the incommensurate spin correlations 
in this material caused by the dominant ferromagnetic 
coupling between the chains. Let us also briefly comment 
on very recent neutron scattering results on powder sam- 
ples of LiCuSb04 [38] . The powder averaging prevents a 
detailed analysis, however, the concentration of spectral 
weight at the incommensurate wave vector qi ^ 0.477r 
points again to a strong frustration a ~ — 2 which is fully 
consistent with an analysis of susceptibility and specific 
heat data [38]. 
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